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Abstract. - We consider fractures in a stratified composite material with solid layers separated 
by thin slices of extremely soft matter. Viscoelastic effects associated with the soft layers are 
taken into account via the simplest model for weakly cross-linked polymers. We find that 
certain small cracks running along layers take a new "trumpet" shape quite different from 
previously known shapes. 



Introduction. - Nacre is a composite material with solid layers (aragonite) separated by 
thin slices of soft organic matter (proteins) (Fig. JlJ). It has a spectacular toughness where 
a fracture propagates normal to the layers This toughness can be explained from 

the absence of stress concentration in this structure J5|,[|. In the present note, we consider a 
different problem, where the fracture plane is parallel to the layers. The fracture properties 
have some similarity with those of a smectic liquid crystal, which have been discussed under 
the name of lenticular fracture Our aim here is to consider the dynamics, when the soft 

layers behave like a viscoelastic fluid. But we shall start with a reminder of the statics, using 
a simple scaling argument. 

The crucial feature is that, for a fracture cavity of size X along the layers, the elastically 
distorted zone has a size Y (perpendicular to the layers) which does not scale like X, but 
rather like 

Y~X 2 /l (1) 
for small cracks (X <C I), where I is a characteristic length 

I 2 = Kb/ Eq. (2) 

Here Kb is a bending modulus and Eq (= eE s ) is an elastic modulus associated with soft 
layers. It is emphasized that Eq. ^) always holds at the scaling level. (For larger cracks 
(X > I) the relation (pi) turns back to X ~ Y; this regime will be discussed elsewhere H). 
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perpendicular fracture 




parallel fracture 



Fig. 1 - Narce-type structure of materials: hard layers (elastic modulus Eh, thickness dh) are glued 
together by soft layers (modulus E s , thickness d a ). Cracks in the y — z plane and the x — z plane are 
called a perpendicular and parallel fractures, respectively (the 2-axis is perpendicular to this page). 



For stratified materials, the following relation has been shown 

l = d/«/I 



(3) 



(4) 



where 

_ Es_ dh_ 

Eh d s 

in the small e limit (see Fig. [l] for notations). In the case of nacre, the small crack condition 
X <C I is rather severe (I ~ 50d). Thus, the following considerations might be more practical 
for some artificially synthesized layered composites where the soft part is a weak gel, Eq is 
very small and / 3> d. 

The small crack condition X <C I makes the situation different from smectic liquid crystals 
where we have Y 3> X because I corresponds to an atomic scale (X/l 1 in Eq. (El)) 
In the present case, for X/l <C 1, we have the anisotropic strain field is distributed 

widely in the x direction compared with in the y direction (see Fig. ^| below) . 

The potential energy (per unit length in the z-direction) of the crack is 



F - X s ( — ) XY- auX + GX 



(5) 



Here, u is the displacement in the y-direction, a the pulling stress (along the y axis) and 
G the fracture energy (energy required to create a new unit area). (It has been shown that 
the dominant component in strain and stress tensors are indeed the y components (it and a, 
respectively) Q). Note here the first elastic term can be equally expressed as Eq{u/Y) 2 due to 
Eq. ([!]) ; this term actually results from a local balance between these two elastic contributions 
andEq. (@) 

originates from this balance condition |7[]. Minimizing F with respect to u, then, 
F(X) has a maximum defining the onset of fracture. At this critical of fracture we have 



X, a ~ X' 



G 



(0) 



where the last equation announces that the product au gives the fracture energy. These 
fractures are very different from conventional parabolic fractures in linear elastic fracture 
mechanics jl0|]: u ~ X 1 / 2 and a ~ X -1 / 2 . In a more precise analysis, we have obtained forms 
for u(x,y) and a(x, y), M which are consistent with Eqs. (0) and ([!]). 
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Fig. 2 - Conceptual (rough) image of strain (stress) and deformation fields. Solid line and broken 
lines are positive and negative deformation contours. Magnitudes of stress or strain are indicated by a 
gray scale (the darker area, the larger stress). Smaller the distance from the tip, the more anisotropic. 
Beyond a distance larger than I, recovers the usual isotropic behavior. (Fracture studied in this paper 
is limited: X < I). 



Eq. (|l]) states that, at the scaling level, a point away from the tip by a distance X along 
the x axis has the same order of strain ( or stress ) with that by a distance Y along the y axis, 
when Y ~ X 2 /I; the situation can be conceptually represented as in Fig. 0. 

We now proceed to the dynamics via a complex modulus of the form W 

fi(<j) = E + (E oo -E )-^- (7) 

1 + iWT 

where the ratio Eoo/Eq = A is assumed to be large as it is in a weakly cross-linked system. 
Here, Eq is related to a small modulus associated with weak cross-links while E^ to a large 
modulus originating from entanglements. Eq. (j^) is a result of taking viscoelastic effects in 
soft layers through a complex modulus of the same form (but with Eq and E^ replaced by 
E s and XE S ) §. 

Crack shape and fracture energy. - We consider a parallel crack propagating with a con- 
stant speed V. When V is smaller than a sound velocity (~ ^Eo/p), the equation of motion 
of the density p and the local velocity v at the scaling level, pDv/Dt = — Vcr, reduces to a 
static equation: Vcr = 0. This is even true in our linear rheological model. In addition, the 
stress components must also satisfy a compatibility equations, which directly result from the 
definitions of strain fields as derivatives of deformation fields. But these geometric conditions, 
again, have the same scaling structures for our linear rheological model. Thus, the scaling re- 
lation for the stress in Eq. (^), a ~ 1/X ~ 1/yT, remains unchanged even in our viscoelastic 
model. These observations will be more precisely addressed elsewhere || . 

Another important observation here is the scaling identification of a distance X along the 
x axis from the tip and the frequency u> via the speed V: 

X ~ V/uj; (8) 

small distances correspond to high frequencies while long distances to low frequencies — the 
farther away from the tip, the more time for relaxation. In addition, from Eq. ffl), the 
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same magnitude of strain with the point specified by the above X is developed at a point 
separated from the tip by a distance Y (~ X 2 /l) along the y axis; these points have the 
same time (1/w ~ V/X) to relax. Thus, a distance Y from the tip sees the same frequency 
frequency u> but a slower propagating speed V y ~ VX/l (Y ~ V y /u>). In other words, we can 
imagine that Fig. ^| correspond to a sequential propagation (from the center to outwards). 
We also note that, seen from a new coordinate (x',y') = (x 2 /l,y), the system returns to an 
isotropic system. For example, Eqs. ( jq) a re changed into the conventional parabolic form: 
cr - 1/VX 1 ~ 1/VY 7 and u ~ 4X 1 ~ VY 7 etc. 

Our viscoelastic model in Eq. (Q) has three regimes depending on frequencies: (I) at 
small frequencies (lot <C 1/A), it is like a soft solid with a small modulus n(u>) ~ Eq, (II) at 
intermediate frequencies (1/A <C lot -C 1), it is like a liquid with viscosity 11(00) ~ iwrj, and 
(III) at high frequencies (1 <C ut), it is like a solid with a large modulus ~ £oc- Due 
to the correspondence between a distance and a frequency in Eq. (@), a fracture can be thus 
spatially divided into three regions (Fig. ||): 

(I) XV t < X: soft solid (modulus £ ) 

(II) Vt < X < AW: liquid (viscosity 77) (9) 

(III) d < X < Vt: hard solid (modulus £Joo) 

/ 

L 

III. SS 

II. L \ 
I. HS~\ 



d 

Vt 
XVx 



Fig. 3 - Three spatial regimes for our viscoelastic model. Due to the scaling relation X 2 ~ Yl, the 
regions appear as anisotropic. The smallest region defined by d comes out for a continuum theory. 
Note that for a crack smaller than in this figure (Vt < L < AW) only the regions I and II are 
developed while for an even smaller crack (d < L < Vt) only the region I is developed; this figure 
corresponds to a fully developed crack (L > AW). Note also that the degree of anisotropy of the 
boundaries separating regions are subject to the ratio of Vt or AW to the length I; the smaller the 
ratio, the more anisotropic. 

We consider a small crack size L with L <§; AW <C I (larger crack sizes (L > !) will 
be discussed elsewhere 0); then, all the three regions (I)- (III) are fully developed. The 
soft-solid region (I) corresponds to low frequencies, and thus, to the static limit; in this 
region (I ^> X S> AW), Eqs. (H) hold. In the liquid zone (II), the stress field scales as 
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a ~ unu/Y ~ nVlu/X 3 and at the same time it should scale as <r ~ Y~ x / 2 ^ \/X even in 
this liquid region as stated above. Thus, the strain should scale as u ~ X 2 and the product 
as era ~ X. The coefficients can be determined by the matching at X ~ Al^r (the boundary 
between I and II): 

GJf 

era ~ (liquid zone). (10) 

XVt 

In the hard-solid region of E x , we find a ~ 1/y/Y and u ~ v 7 ! 7 as in Eqs. @ but with 
era ~ Go (via the same manner as in deriving Eqs. @). Here, Go is associated with the 
hard solid appearing near the tip. Matching this latter product era with that in Eq. (|l(]) at 
X ~ Vt, we find G ~ AGo- The overall separation energy G for a fully developed crack is 
enhanced from Go associated with local precesses near the tip. Note here that this expression 
for G is valid for d < Vt in our continuum theory. 

The crack shape resulting from this analysis is summarized as follows: 



X for XVt < X 

X 2 for Vt < X < XVt 

X for X < Vt 



(11) 



It is just like a trumpet with a lenticular edge (Fig. Q), as has been suggested by the name of 
the model, but different from previously known shapes; it is not similar to the conventional 
parabolic form nor an isotropic parabolic trumpet predicted |pd| and observed fl2|| in certain 
polymer systems (Fig. ^). 



(a) lenticular trumpet 
X 



(b) parabolic trumpet 





Fig. 4 - Lenticular and parabolic viscoelastic trumpets. 

We complete our arguments by considering smaller fractures. When Vt < L < XVt, only 
the hard-solid and liquid region are present; the soft solid has yet to develop. In this situation, 
the fracture energy is given by Eq. ( [To] ) at X = L: GqL/ (Vt); the toughness decreases with 
velocity. When L < Vt, only the hard-solid region is developed and the fracture energy is 
given by Go- Thus, with increase in V, the fracture energy starts from a larger plateau value 
AGq, and then decreases to reach a smaller plateau value Gq: 



f AGo 

G(V)~{ G L/{Vt) 
{ Go 



for d/r <V< L/(Xt) 
for L/(Xt) <V < L/t 
for L/t <V 



(12) 
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This behavior can be confirmed more precisely from a general formula: 

G(V) 



Go 



E 





1 


/^Im 




/ u 



(13) 



which can be analytically calculated for the present model (just as in the previously known 
isotropic case fl3|]). It should be emphasized here that this formula is unaltered even in our 
anisotropic materials, which is by no means trivial; Eq. jl3| ) here can be shown in the following 
manner. We start from a relation: 

G(V)V ~ J dx J dyoe ~ J dXYoxex- (14) 

In order to estimate ex we again use the dimensional identification in Eq. (ih: ex — e^ ~ 
uj<j u / h(uj) ~ luox/ and 

G(V) ~ / dcj^- ~ ^Gq / ^- (15) 

Here, we have used a more precise form of Eq. (^): erx ~ \/E~^Go~/Y ||. Since the real and 
imaginary part of l/fj,(u>) are even and odd functions, respectively, we arrive at Eq. (|l3|). 

Conclusion. - In this paper, we present a physical picture for fractures in nacre-type 
materials via scaling arguments. Viscoelastic effects for parallel fractures are taken into ac- 
count via the simplest viscoelastic model for weakly cross-linked polymer. We expect that for 
slow crack-propagation speeds (I ^> XVt) a small crack (L -C I) takes a new trumpet shape 
different from previously reported shapes; in the opposite limit will be discussed in a separate 
paper || . The overall fracture energy G is found to decrease from a larger plateau value XGq 
to a smaller plateau value Go with increase in velocity V where Go is associated with local 
precesses near the tip. 
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